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Abstract 



In numerical algebraic geometry witness sets are numerical representations of positive 
dimensional solution sets of polynomial systems. Considering the asymptotics of witness sets 
we propose certificates for algebraic curves. These certificates are the leading terms of a 
Puiseux series expansion of the curve starting at infinity. The vector of powers of the first 
term in the series is a tropism. For proper algebraic curves, we relate the computation of 
/--^ • tropisms to the calculation of mixed volumes. With this relationship, the computation of 

tropisms and Puiseux series expansions could be used as a preprocessing stage prior to a more 
expensive witness set computation. Systems with few monomials have fewer isolated solutions 
and fewer data are needed to represent their positive dimensional solution sets. 
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\ 1 Introduction 

CN 

Solving polynomial systems numerically used to be restricted to finding approximations to all 
isolated solutions. Via numerical algebraic geometry [63} [M] we are now able to provide numer- 
ical representations for all solutions and for all dimensions. The development of the methods of 
numerical algebraic geometry coincided with upgrades [61] to the software PHCpack [75] . Its 
blackbox solver computes isolated solutions via polyhedral methods. The aim of this paper is to 
bring polyhedral methods to the foreground and to show how they may lead to finer representa- 
tions of solution sets. 

We consider as given a polynomial system with at least as many equations as unknowns, and 
we want to determine whether a proper algebraic curve is a solution. An algebraic curve is proper 
if it is not contained in any higher dimensional solution set. Moreover, we will assume that the 
proper algebraic curves we are looking for are of multiplicity one. What we want to determine of a 
solution curve is a certificate of its existence. The certificate should be small and straightforward 
to use in computer algebra systems. 
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For the solution to this problem we propose polyhedral methods to find the leading terms 
of a Puiseux series expansion [54J of the solution curve. Verifying the existence of a curve as a 
solution of the given polynomial system requires then only the formal substitution of the Puiseux 
series into the given system, a routine task for computer algebra systems. The main source of 
inspiration for this approach was found in the emerging field of tropical algebraic geometry [33] , 
[55], [67] • To prove that a polynomial system from celestial mechanics has only finitely many 
isolated solutions, polyhedral methods were applied in [30]. This paper could be seen as an 
attempt to make such proofs automatic. 

Working with series expansions is a hybrid form of computation, combining symbolic and 
numerical approaches. Some algebraic curves may actually have finite series expansions and 
in the extreme case even consist of only one leading term. When considering systems with 
natural parameters, such solution curves may be more useful than other solutions and it may be 
worthwhile to look for such solutions first. Our goal is to develop a polyhedral method that will 
focus on computing one dimensional solutions. 

After defining witness sets and outlining our problem statement, we will introduce tropisms 
for binomial systems. Proper algebraic curves defined by binomial systems admit a very explicit 
solution. To examine the normalization of tropisms we consider deformations of witness sets. The 
link between tropisms and mixed volumes passes through the second theorem of Bernshtem [Sj- 
We show how to compute tropisms to proper algebraic curves defined by square systems via a 
special lifting. 

In the spirit of the theme of the conference, held in honor of Andrew Sommese, 22-24 May 
2008, this paper is on the interactions of classical algebraic geometry — in particular the Puiseux 
series, the roots of polyhedral methods and tropical algebraic geometry — and the new field of 
numerical algebraic geometry. Detailed descriptions of algorithms and their implementations are 
still a work in progress. 

Acknowledgements. Preliminary versions of the ideas in this paper were presented at four 
conferences held respectively in Tobago, Notre Dame, Hong Kong, and Vancouver. The author 
thanks all organizers of these meetings for their invitations and the participants for their feedback. 

2 Witness Sets and Newton Polytopes 

In numerical algebraic geometry, homotopy continuation methods ( [40] . |48| . |64j ) manipulate 
algebraic sets. The key data representation for algebraic sets is a witness set (first defined as 
witness point set in [60], see also [64])) defined below. 

Definition 2.1 Given a system f(x) = 0, we represent a component of f _1 (0) of dimension k 
and degree d by a witness set which consists of f , k general hyperplanes L and d generic points 
in f~ 1 (0)nL. 

A generic choice of coefficients for the k hyperplanes implies that all solutions will be isolated 
and regular, unless the component has a multiplicity higher than one. Definition 12.11 incorporates 
the theorem of Bertini, see |64] for its application in numerical algebraic geometry. Bertini [6] is 
a more recent software system for numerical algebraic geometry. A dictionary compares in [62] 
witness sets to lifting fibers [27] . [28] in a geometric resolution [26]. 
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We call a system square if it has as many equations as unknowns. After adding linear equations 
to a square system, we |59j add slack variables in an embedding to make the overdetermined 
system square again. Using a flag of linear spaces, defined by an decreasing sequence of subsets 
of the k general hyperplanes, 

L = L k D L fc _x D • • • D Li D L = 0, (1) 

we move solutions with nonzero slack values to generic points on lower dimensional components, 
using a cascade of homotopies [59]. By default, the top dimension k in the cascade starts at n — 1, 
but if we know the system has only solution curves we start with k = 1. According to [25] . the 
dimension of an algebraic variety can be determined in polynomial time. 

The cost of the homotopy cascade algorithm is determined by the number of paths, starting 
at the solutions of the top dimensional system. By the embedding, this top dimensional system 
has only isolated roots and thus one may apply any solver to compute those isolated roots. To 
count the number of solutions for square systems with generic coefficients we look at the Newton 
polytopes. The Newton polytopes are spanned by the exponents of the monomials which occur 
with nonzero coefficients in the system. We formally describe sparse polynomials as follows. 

Definition 2.2 Given a polynomial 

/(x) = ^c a x a , c a ^0, x a = xf4 2 (2) 

a<=A 

the support A of / collects only those exponent vectors with nonzero coefficient in /. The convex 
hull of A is the Newton poly tope of /. We denote the inner product by (•, •). To define faces d v P 
of P we use a support function p: 

p(v) = min (x, v) so d v P = { x £ P | (x, v) = p(y) }. (3) 

The equation (x, v) = p(v) determines a supporting hyperplane for the face d v P. A vector v 
perpendicular to a fc-dimensional face d v P of P lies in an (n — A;)-dimensional cone. 

In defining support functions as in ([3]), we choose the minimum instead of the maximum 
convention, opting for inner rather than outer normals. These minimum and maximum conven- 
tions correspond to letting a deformation parameter respectively go to zero or to infinity. In [24] , 
Newton polytopes arise as compactifications of amoebas [35], obtained by taking logarithms of 
the variety [7.]. Algorithms to compute amoebas are presented in [68] . We refer to [69] and [76] 
for references on polytopes. For a system f(x) = 0, we collect the Newton polytopes of the tuple 
of polynomials f = (/i, f 2 , ■ ■ ■ , /„) in the tuple P = (P 1 , P 2 , . . . , P n ). 

Definition 2.3 Given a tuple of Newton polytopes P = (Pi, P%, . . . , P n ) we define the mixed 
volume V^(P) via the formula 

V n (P 1 ,P 2 ,...,P n ) = Y, Pl(v) K-l(3 v P 2 ,...,5 v P n ), (4) 

gcd(v) = 1 

where p\ is the support function for P\, Vectors v are normalized so the components of v have 
their greatest common divisor equal to 1. 
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Because polytopes are spanned by only finitely many points, only finitely many vectors v will 
yield a nonzero contribution to In the special case: P = P\ = Pi = ■ ■ ■ = P n : V^(P) = 
n ! volume (P). Handling the mixed case, when several but not all polytopes are repeated, is 
important for efficient algorithms, developed in [21], [51], and |47j . Mixed volumes are at the core 
of classical geometry [58J and have applications to tomography |23j . 

Theorem 2.4 (Bernshtein Theorem A [8]) The number of roots of a generic system equals 
the mixed volume of its Newton polytopes. For any system, the mixed volume bounds the number 
of isolated solutions in (C*) n , C* = C \ {0}. 

In the same paper, Bernshtein gave in his second theorem precise conditions for the mixed 
volume to be sharp. Because of these conditions, we speak of count [12] instead of bound, see 
also |56l I57j. We will formulate the second theorem later, but now we have introduced enough 
terminology to formulate our problem statement. The proof for Theorem 12.41 given in [8] is 
constructive and served as basis for polyhedral method in [71] to find all isolated solutions. A 
more general polyhedral homotopy homotopy method was developed in [31] . see also [66J. The 
authors of [36] address the complexity of Bernshtein's first theorem. The development of recent 
software for polyhedral homotopies is described in [29] and [38] . 

We list three problems with the current use of witness sets. First and foremost, as we add 
hyperplanes to the system, the polyhedral root count increases for sparse systems. For example, 
for a benchmark problem like the cyclic 8-roots system, adding one hyperplane raises the mixed 
volume from 2,560 to 4,176. For cyclic 12-roots, the same operation brings the mixed volume 
from 500,352 to 983,952. While the cascade will give us all solutions at the end, if we are only 
interested in the curves, we do not want to compute start solutions to the isolated roots. The 
second problem concerns symmetry. Many polynomial systems have obvious permutation sym- 
metries and with polyhedral methods we can setup symmetric homotopies |72j to compute only 
the generators of each orbit of isolated roots. However, extending the symmetric polyhedral ho- 
motopies of [72] to deal with positive dimensional solution sets conflicts with the current witness 
set data representation. Thirdly, users of numerical algebraic geometry methods require and need 
guarantees for the results to be correct. Once we have a witness set, any path tracker may be 
used to sample points on the solution set. While path trackers are standard in numerical anal- 
ysis, users unfamiliar or uncomfortable with floating-point computations require exact answers. 
Additionally, to determine the degree of a solution set correctly, we must be able to certify that 
all solutions have been found. 

3 Proper Algebraic Curves defined by Binomial Systems 

A binomial system has exactly two monomials with a nonzero coefficient in every equation. We 
consider n — 1 equations in n variables. By limiting the number of equations and restricting 
to nontrivial solutions in (C*) n we are reducing the complexity of the problem. However, even 
already for binomial systems, the complexity of counting all isolated solutions is #P-complete |13j . 

For example (n = 3): 
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The system at the right of ([5]) is a normal form of the system. For general n, we can always write 



a binomial system as a tuple of equations of the form x s 



c. 



Writing a system in this normal 



form removes trivial solutions with zero coordinates. For a binomial system in its normal form, 
we collect all exponents in a matrix A G 1 ) xn and its coefficients in c G (C*)™ -1 . Continuing 
the example, we have: 



.4 



-1 
1 4 



rank(A) = 2, v 



: Aw = 0. 



(6) 



The vector v in the kernel of A will determine the shape of the solution curve. 

For A: rank(A) = n — 1, there is a unique vector v in the kernel. A binomial system of 
n — 1 equations in n variables will have have a proper (i.e.: not contained in any other higher 
dimensional solution set) solution curve in C* if and only if the rank of the exponent matrix A 
is n — 1. Although a normal form of a binomial system is not unique, the vector v for a proper 
algebraic curve is independent of the choice of a particular normal form. We will show that proper 
algebraic curves have solutions of the type Xk = c\ t t Vk , with ct G C*, k = 1, 2, . . . , n. 

To simplify the system, we use v to define a unimodular matrix M (det(M) = 1) and a 
coordinate transformation (called a power transformation in [11]) denoted by x = y M : 



M 



+4 1 
-4 1 
-10 



AM 



XI 


= vl 










X 2 


= Vi 








1 


x 3 


= Vi 



(7) 



After applying the coordinate transformation, defined by M, we can divide out the variable y\ 
and we obtain a system of two equations in two unknowns: 1/2 and 2/3. For our example, we find 
one solution 



(2/2 



T2-,V3 



-6) to 



y 2 1 1 = 2 



112 



(8) 



We rename the free variable y\ to t and use the values found for 1/2 and 7/3 in the representation 
of the solution curve, see the left of ([9]). To compute the degree of the curve, we add a random 
hyperplane and substitute the expression found for the solution: 



-6i+ 4 

12 b 



-1 



,-1 



X 2 Xl 



2 
-3 



(9) 



70 + 71^1 + 72^2 + 73^3 = 



The coefficients 70, 71, 72, and 73 are random complex numbers. After substitution the left of © 
into the right of Q we obtain 70 + 71 (— 6i +4 ) + 72(— j2~t~ A ) + 73^ _1 - Clearing denominators, we 
find a polynomial in t of degree 8. 

Given the leading term of the Puiseux series expansion for an algebraic curve, the degree of 
the curve follows from the leading exponents of the Puiseux series. We can formalize this in the 
following proposition: 
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Proposition 3.1 Consider a proper algebraic curve defined by x = c, with v: Av = 0. Let M 

be the unimodular matrix in the transformation x = y M that eliminates y\ . Let B be the matrix 
obtained from removing the first zero column of AM . Then the degree of the curve equals 

n n 

| det(£?) x | max Vi — min V{ \ . (10) 

i=l i=l 

Proof. The form of the solution curve is = ot^t Vh where the coefficients are the nonzero roots 
of y AM = c. The number of roots equals | det(B)|. The curve will have as many components 
as | det(P)|. Because every component has the same degrees in i, every component will have the 
same degree, so it suffices to compute the degree of one component. 

To compute the degree of one component, we reduce it to as many isolated points as its degree, 
intersecting it with a hyperplane with random coefficients. Substituting the form of the solution 
component into that hyperplane yields a polynomial in t. To clear denominators we multiply by 
the most negative exponent or we divide out trivial solutions by subtracting the lowest positive 
exponent of v. □ 

We point out that the formula (llOj) is not invariant to unimodular transformations. Consider 
for example the plane curve defined by f(x\,X2) = x\X2 — 1 = with v = (+1,-1). The 
unimodular coordinate transformation defined by v: x\ = y\, Xi = y^ X y2 reduces / to y% — 1 = 0. 
In the new ^-coordinates, we now have a line (yi, 1). 

For a proper algebraic curve defined by a binomial system we can provide an exact certificate 
for its degree, independently of the choice of the coefficients. The vector v we computed in the 
example above is an example of a tropism. We define tropisms for general polynomial systems in 
the next section. 



4 Tropisms and Initial Forms 

Puiseux series occupy a central role in the study of algebraic curves. The leading exponents of 
the series are called tropisms, defined by the Newton polytopes of a polynomial system. 

Definition 4.1 Consider f(x) = with Newton polytopes in (Pi, P2, . . . , Pjv). A tropism is a 
vector perpendicular to one edge of each Pj, for i = 1, 2, . . . , N. 

Our definition of tropisms differs from the usual one in the literature (|39j. |43j ) where all 
coordinates in the tropism (or critical tropism in singularity theory [39]) are required to be 
positive. This requirement is natural if one looks for one point on a solution curve in affine space 
and then shifts that point to the origin. The only normalization we will require is on the sign 
of the first coordinate. We will provide arguments for this normalization when we consider the 
asymptotics of witness sets in the next section. The other difference in the definition is then 
that we may have tropisms pointing out isolated solutions at infinity, solutions that do not give 
rise to an initial term in a Puiseux series expansion. Using a more refined terminology we could 
give the v in Definition 14.21 the name pretropism, as a part of a tropical prevariety. A tropical 
prevariety [10] corresponds to the intersection of the normal cones of the polytopes. 
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The edges perpendicular to a tropism are Newton polytopes of an initial form system which 
may have solutions in (C*) ra . 

Definition 4.2 Let v£Z"\ {0} and /(x) = V c a x a . Denoting the inner product by (•,•), the 

aeA 

initial form of f in the direction v is 

in v /(x) = c a x a with m = min{ (a, v) | a € A }. (11) 

a. £ A 
(a, v) = m 

Let f = (/i,/2, • • • >/iv) be a tuple of polynomials. For v E Z n \ {0}, the initial form system 
in v f (x) = is denned by the tuple in v f = (in v /i, in v /2, • • • , in v /iv)- 

Although we define initial form systems via tropisms, vectors perpendicular to edges of the 
Newton polytopes, the number of points in the support of each initial form is at least — but not 
exactly — two. Although the initial form system is thus not necessarily a binomial system, we 
can always eliminate one variable. 

In [11] and [37], systems supported on faces of Newton polytopes are called truncated systems. 
This terminology refers to the process of substitution the power series and then selecting those 
terms that correspond to the lowest power in the variable of the series. We prefer to use initial 
forms because of the relationship with term orders widely used for Grobner bases [65], |69j . 

Of special importance are the conditions for which the mixed volume is sharp, formulated in 
the tropical language. 

Theorem 4.3 (Bernshteln Theorem B |8j) Consider f(x) = 0, f = (/i, f%, . . . , f n ), x = 
(xi,X2, • • • ,x n ). If for all tropisms v; in v f(x) = has no solutions in (C*) n , then f(x) = has 
exactly as many isolated solutions in (C*) n as V^(P). 

If there is a positive dimensional solution set, then this set stretches out to infinity and the 
system /(x) = must have solutions at infinity. Of particular interest is the following consequence 
of Bernshtem's second theorem. 

Corollary 4.4 If there are no tropisms v for which in v f(x) = has roots in (C*) n then f(x) = 
has no solutions at infinity. 

The tropisms and the roots of the corresponding initial form systems give the leading term 
of the Puiseux series expansion of the solution curve. The second term in the Puiseux series 
expansion will provide the certificate for the existence of the solution curve. Also the coefficients 
of the second term will be roots of a truncated polynomial system. Once we have those roots, we 
can further grow the Puiseux series expansion symbolically, or apply numerical predictor-corrector 
methods to sample points along the solution curve. In Figure H] (slightly adapted from [1]) we 
sketch the idea for computing this certificate. 

The computation of the second term in the Puiseux series goes along the so-called Newton- 
Puiseux method, as outlined in the proof of the theorem of Puiseux [75], see also [UJ. In [2] 
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Newton polytopes 

I 

compute tropisms O 



solve initial forms O 



compute 2nd term O 



series 



no tropism 
=>■ no root at oo 

no root at oo 
=>• no series 

no series 

no curve 



Figure 1: Computing a certificate for a proper algebraic curve. 



algorithms for Puiseux series for space curves are described and an implementation in CoCoA is 
mentioned. General fractional power series solutions are described in [44] . See [35], [36] and [53] 
for recent symbolic algorithms, and [51], [52] for a symbolic- numeric approach. 

Algebraic curves defined by binomial systems can be solved by one tropism but it may happen 
that one tropism solves a more general system. Consider for example the cyclic 4-roots problem. 



f (x) = ( 



X\ + x 2 + x 3 + x 4 = 
x\x 2 + x 2 x 3 + x 3 x 4 + X4X1 = 
X1X2X3 + X2X3X4 + X3X4X1 + X4X1X2 = 
X1X2X3X4 — 1 = 



(12) 



X1X2X3X4 — 1 = 

There is one tropism v = (+1, —1, +1, —1) which leads to the initial form system in v /(z) = 0: 



in v f (x) 



vt 1 



(13) 



X2 + X4 = [ Xi 

xix 2 + x 2 x 3 + X3X4 + X4X1 = x 2 
X2X3X4 + X4X1X2 = I X3 

X1X2X3X4 — 1 = \ Xi 

The system in v f (y) = has two solutions, and we find the two solution curves: (t, — t -1 , —t, t^ 1 ) 
and (t, — — - 

Note that v = (—1, +1, —1, +1) is a tropism as well for the solution curves of cyclic 4-roots, 
but considering this tropism corresponds to setting x\ = t^ 1 or moving the curve to infinity 
instead of to zero as t goes to zero. We will examine this in greater generality in the next section. 



5 Asymptotics of Witness Sets 

One way to compute tropisms would be to start from a witness set for an algebraic curve in 
n-space given by d points on a general hyperplane cq + C1X1 + C2X 2 + • • • + c n x n = and satisfying 
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a system f(x) = 0. We then deform a witness set for a curve in two stages: 



1. The first homotopy moves to a hyperplane in special position: 

h(xt) = { f(x)=0 

1 (cq + ci^i + • • • + c n x n )t + (cq + cixi)(l — t) = 0, for i from 1 to 0. 



2. After renaming Co + c\Xi = into xi = 7, we let x\ go to zero with the following homotopy: 

h ^ X ' = { xi - jt = 0, for t from 1 to 0. ^ 15 ' ) 



The two homotopies need further study. In the first homotopy (|14p some paths will diverge, 
consider for example f(x\, X2) = x\X2 — 1. Even all paths may diverge if the solution curve lies in 
some hyperplane perpendicular to the first coordinate axis x\ = c with c different from —cq/c\. 

Lemma 5.1 All solutions at the end of the homotopy h(x, t) = of (|14p lie on the curve defined 
by f(x) = and in the hyperplane x\ = —cq/c\. 



Proof. We claim that we find the same solutions to h(x, t = 0) = either by using the homotopy 
in (|14p or by solving h(x, 0) = directly. This claim follows from cheater's homotopy |42j or the 
more general coefficient-parameter polynomial continuation [49]. □ 

The second claim we make in the main theorem below is that we recover all data lost with 
tropisms. For simplest example of the hyperbola X1X2 — 1 = 0: its solution is (x\ = t, x<i = i -1 ) 
and the tropism is v = (1,-1). The lemma below extends the normal form for the Puiseux series 
expansion for plane curves (as used in |75| ) to general space curves. 

Lemma 5.2 As t — > in the homotopy (|15p . the leading powers of the Puiseux series expansions 
are the components of a tropism. In particular, the expansions have the form 

{ xl = c k t v *(l + 0{t)), k = 2,...,n. (16) 

Proof. Following Bernshtem's second theorem, a solution at infinity is a solution in (C*) n of an 
initial form system. For a solution to have values in (C*) n , all equations in that system need to 
have at least two monomials. So the system is an initial form system defined by a tropism. To 
arrive at the form of (|16p for the solution defined by the homotopy (|15p we rescale the parameter t 
so we may replace x\ = yt by x% = t. □ 



Also in the second homotopy, solution paths are most likely to diverge and the directions of 
the diverging paths are defined by the tropisms. The numerical computation of those directions 
can be done by endgames using power series as in [50] and [32] . 
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Definition 5.3 Given a system f(x) = which defines a proper algebraic curve. Consider a 
Puiseux series expansion of the form 



{ 



X\ = t 

x k = c k t Vk + d k t Wk + • • • c k ,d k eC*,v k ,w k eQ,k = 2,... 



n. 



(17) 



Then a certificate for the solution curve consist of the exponents (i>2, ■ ■ ■ , v n ), (1V2, ■ ■ ■ , w n ) £ Q n 1 
and the coefficients (c2, . . . ,c n ), (^2, . . . ,d n ) G (C*) n_1 . 

The tropism v (or pretropism) shows there are solutions at infinity, but solutions at infinity 
could be isolated. In that case the tropism v is still a certificate (but then more like a death 
certificate) for the lack of sharpness of the mixed volume to count all isolated roots. The expo- 
nents w of the second term in the series show the solution is part of a curve. Arguing in favor of 
extending the term tropism rather than reserving it only for the leading exponents in the series, 
we point out that a certificate consists of a pair of vectors: (v, w) and a corresponding pair of 
solutions (c, d), respectively of an initial form system and a truncated system. As solutions of 
overdetermined polynomial systems, the coefficients (c, d) can be certified [18] by a-theory [9]. 

Relating the data (v, w) and (c, d) to a witness set, we note that the data corresponds 
geometrically to cutting the curve with a special hyperplane x\ = 0. There is the risk of missing 
curves in coordinate hyperplane x\ = c, for c G C* and there may also be singular solutions 
occurring for x\ = 0. Except for these two difficulties, the certificate provides a predictor to 
sample the solution curve. The capability to sample a solution curve is intrinsic in the definition 
of a witness set, but to verify this capability one needs to compute at least one step of Newton's 
method at at least one point of the witness set. Substituting the certificate in the original system 
— with and without the second term and watching the degree in t increase — is a more elementary 
operation. 

On the one hand, we can view the certificate as a very special witness set, obtained by 
intersecting the curve with the hyperplane x\ = 0. On the other hand, we can view the certificate 
as a very special lifting fiber, where one free variable x\ is specialized to zero. 

Lemma 5.4 The certificate of (I17p may be written in the form 



where a k ,Pk £ C* and v k , G 2 and v\ is the smallest natural number to clear the denominators 
in the series for the other components Xk, k > 1. Then the degree of the branch is determined by 



where R is the set of initial roots of the initial form system in^f (x = y ) = 0, forM a unimodular 
matrix with first column equal to v = (y\, U2, ■ ■ ■ , Vn)- 

Proof. The existence of v\ follows from the definition of Puiseux expansions. Suppose there would 
not be a smallest v\ to clear denominators, then we could make a plane curve for which we could 
not clear denominators. Proposition 13.11 is generalized in two ways: 




(18) 



fj^R x | max Vi — min vi 



1=1 4 = 1 



(19) 
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1 . The initial form system in^f (y ) = is no longer a binomial system for which we can count 
the number of roots via a determinant. Instead we now solve m„f(y) = and collect the 
roots in R. 



2. We have to show that only the leading terms in the Puiseux series expansion determine the 
degree of the solution branch. Consider t — > and consider the solutions in the hyperplane 
x\ = t Ul . By the form of the expansion we cannot have more solutions for t > then we 
would have at t = 0. At t = 0, only the leading terms matter in defining the initial form 



The key point of the two homotopies we considered above is the argument for the normalization 
of the first coordinate v\ of the tropisms to v i > 0. Moving the degree many points on the curve 
on the hyperplane x\ = 7 to x\ = with x = t is equivalent to moving those points to infinity 
with x = t . By equivalent we mean that we do not obtain any new information about the curve 
by considering also tropisms with first coordinate of the opposite sign. 

Theorem 5.5 Given a proper algebraic curve defined by a polynomial system f(x) = ofn equa- 
tions in n unknowns. Assume that the curve does not lie in a coordinate hyperplane perpendicular 
to the first coordinate axis. Then the degree of the curve corresponds to the number of certificates 
times the degree of each certificate. 

Proof. If we consider the two asymptotic homotopies in sequence, then we may lose solutions in 
the first move, by putting the cutting hyperplane perpendicular to the first coordinate axis, while 
in the tropisms found in the second stage lead then to higher degrees. The claim of the theorem 
is to show that a potential loss of witness points is made up by the increase in the degree of the 
Puiseux series at the very end. 

Executing the asymptotic homotopies in sequence is equivalent to combine them into one 
homotopy: 



Observe that the Newton polytopes of f are not deformed by this homotopy. So the structure of 
the space at infinity remains invariant as well. Applying polyhedral endgames [32J (see also [50J), 
the witness points will end at roots of initial form systems. So every witness point corresponds 
to one initial root and we cannot have more witness points than we have initial roots. 

Using a weighted projective space |14j . |15j . |71j . determined by the tropisms, we can extend 
the initial roots to roots for t > 0. So the initial roots give rise to witness points. □ 

We end this section with a note on the complexity of the certificates. In the best case, the 
initial form systems are binomial, while in the worst case the number of monomials in the initial 
form systems is of the same order of magnitude as the original system. But even in this worst 
case, the number of variables drops by one and this drop may be enough to get a more tractable 
problem. 



system in iy f(y) = 0. 



□ 




for t from 1 to 0. 



(20) 



11 



6 Tropisms and Mixed Volumes 



Degenerating witness sets is an effective but not an efficient way to compute tropisms. Relating 
to mixed volumes, we indicate how to compute the tropisms directly from the Newton polytopes. 
Before we explain the lifting algorithm to compute mixed volumes we point out that Gfan [M] 
using the algorithms of [10] is more appropriate to compute tropisms. Our point of relating 
tropisms to mixed volumes is to investigate the connection between the generic number of isolated 
roots and the degrees of the solution curves of sparse polynomial systems. 

Let A = (Ax, A2, . . . , An) be the supports of f(x) = 0. Following [19], [31], we summarize the 
mixed volume computation in three stages: 

1. Lift a G A using a function to: u>(Ai) C R n+1 . 

n 

2. The facets on the lower hull of the Minkowski sum ' s ^^uj{Ai) spanned by one edge of each 



of uj(Ai) define mixed cells C . 
3. The mixed volume is V n (A) = Vol(C). 

C is mixed 



i=l 



By duality [73], mixed cells are defined by inner normals perpendicular to edges of the polytopes. 
These inner normals are tropisms with positive last coordinate. 

Polyhedral homotopies [31] follow the computation of the mixed cells and the lifting function 
defines the powers of the new parameter t. In particular, polynomials are lifted as 

/(x) = ^c a x a - /(x,t) = ^c a x a ^(a) (21) 

aeA aeA 

using the same lifting function to as before. 

Now we look for solution curves of the form (|16j) . We observe that since x\ = t, as lifting 
function we define t as r(a) = deg(x a ,xi) = a\ as lifting function. The lifting r will work for 
systems of n — 1 equations in n variables. However, for systems with n equations, if we use t = x±, 
then we will have too few variables. Therefore, we use a slack variable z in the lifting, and define 
r as 

r : C[x] -> C[z,x] 

Y a _ «i «2 a n ._. _/' Y a\ _ r a 2 a n j.a-i 

-A. — ju -y jl> 2 n V / — 2 n 

where r is some random exponent. For exponent vectors a, we define r as r(a\,a2, ■ ■ ■ ,a n ) = 
(r, a,2, ■ ■ ■ , OLni a l)- Note that r = w(a), the usual random lifting used to calculate mixed volumes. 
If we take t = X\, z = t, r = uj(a), then we get i w (a)x2 2 • • • x^x^ 1 = x a t ul (a), which is the random 
lifting commonly used in polyhedral homotopies [31] . With the slack variable we have again as 
many variables as equations and we can apply our mixed volume calculators. 

We claim that the tropisms to proper algebraic curves are in one-to-one correspondence with 
those inner normals to the mixed cells of the mixed subdivision induced by the lifting u for which 
the z-component is zero. 
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Proposition 6.1 Let S be the set of mixed cells for an n-tuple A in n-space lifted with t, intro- 
ducing a slack variable z. Then v = (v z , v x ) is a tropism if and only ifv z = 0. 



Proof. Denoting the inner normals to the cells in S as v = (v 2 ,v x ), points (a, b) in each cell 
satisfy 

r a y z + a x v x = r b v z + b x v x , a = (r a , a x ), b = (r b , b x ). (23) 
We have to show two things: 

1. Each inner normal v with v z = is a tropism. 

If v z = 0, then a x v x = b x v x holds and the tropism is v x because (a x , b x ) span an edge of 
one of the Newton polytopes. 

2. For every tropism we must have that v 2 = 0. 

Given a tropism v x , we have a x v x = b x v x . Rewriting (I23p to solve for v z leads to 

v 2 (r a - r b ) = b x v x - a x v x = 0. (24) 
Since r a and r b are random numbers, v z = 0. □ 

Although this is already a more efficient method than applying polyhedral endgames in the 
asymptotic homotopies on witness sets, for practical purposes one would include the constraint 
v z = already immediately in all feasibility tests the mixed volume calculator does. Ultimately, 
unlike the lift-and-prune approach [19] for mixed volume computation, the complexity of the 
problem of computing tropisms is governed by the shape of the polytopes and the relative position 
of the polytopes with respect to each other. Tropisms lie in the common refinement of cones of 
inner normals to faces of the Newton polytopes and algorithms of [10] as implemented in Gfan [M] 
are recommended. 

Proposition 16.11 relates the number of isolated solutions to components of solutions. Consider 
for example the following system: 



f(x) 



X1X3 — X\ — X3 + 1 = 

x 2 - 1 = (25) 
%3 ~ 7 = 7 G C*. 



The mixed volume of f(x) = equals one. For all nonzero choices of 7 (except for 7 = 1), the 
system will have the isolated solution (1, 1,7). For 7 = 1, we have the tropism v = (1,0,0) and 
the general solution line (xi = t, x-i = l,Xs = 1). A witness set representation for this system 
would need to store two solutions: one solution for the witness set of the line in case 7 = 1 and 
one solution for the path leading to the isolated root in case 7 / 1. 



7 Preliminary Computational Experiments 

The cyclic re-roots problem is a widely known benchmark for polynomial system solvers, see 
e.g. [16] . [20] and [41 j . For those re which are divisible by a square (e.g.: re = 4, 8,9,12), the 
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system is known to have positive dimensional solution sets [3]. The general recipe to formulate 
the polynomial equations for any n in this family is obvious from the cyclic 4-roots system, 
given above in (|12ft . The permutation symmetry in cyclic n-roots is generated by two elements 
(xi,x 2 ,x 3 , ...,x n ) -» (x 2 ,x 3 , . . . ,x n ,xx) and (xi,x 2 , ■ . . ,x n _i,x n ) -> (x n ,xi,x 2 , ■ ■ .,s n -i)- 

The numerical computations reported below were done with PHCpack [70], using a modified 
lifting in MixedVol [22]. An alternative mixed volume calculator is DEMiCs [46J. For the symbolic 
manipulations, Maple 11 was used. We emphasize that the computations reported below are 
preliminary, mainly to illustrate the concepts. 

7.1 cyclic 8-roots 

The program computed 29 tropisms, listed below in Table [Hand Table [2j 
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Table 1: First list of 13 tropisms for cyclic 8-roots, separated in 6 orbits. 
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Table 2: Second list of 16 tropisms for cyclic 8-roots, separated in 8 orbits. 

The last tropism from Table [2j v = (1, 0, —1, 1,0, —1, 0, 0), the initial form system is 

x 3 + x 6 = 
x 2 x 3 + x 5 x 6 + X 6 X 7 = 
X5X6X7 + X6X7X8 = 
£3X4X5X6 + X 5 X 6 X7X 8 = 
£2X3X4X5X6 + X 3 X4X 5 XqX 7 = 
X2X3X4X5X6X7 + X3X4X5X6X7X8 + XqXjXsXiX 2 X 3 = 
X2X3X4X5X6X7X8 + X 5 XsX 7 X$XlX 2 X 3 = 
XiX 2 X 3 X<iX5XqX 7 X8 — 1 = 



in v f (x) 



(26) 
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To reduce the system into a simpler form, we perform the coordinate transformation, dividing 
out yi in the initial form system: 



Xi 
X2 



yi 

2/l2/2 

£3 = 2/^2/3 
xa = y\yi 

xg = yi l y§ 

xs = 2/12/8 



in v f(y) 



2/3 + ye = 

2/22/3 + 2/52/6 + 2/62/7 = 

2/52/62/7 + 2/62/72/8 = 

2/32/42/52/6 + 2/52/62/72/8 = 

2/22/32/42/52/6 + 2/32/42/52/62/7 = 

2/22/32/42/52/62/7 + 2/32/42/52/62/72/8 + 2/62/72/82/22/3 = 

2/22/32/42/52/62/72/8 + 2/52/62/72/82/22/3 = 

2/22/32/42/52/62/72/8 - 1 = 



(27) 



To solve the overconstrained initial form, we introduce a slack variable s and generate eight 
random numbers 7^ £ C*, k = 1, 2, . . . , 8 to multiply s with. Then we solve 

2/3 + 2/6 + 7i« = 
2/22/3 + 2/52/6 + 2/62/7 + 72 s = 
2/52/62/7 + 2/62/72/8 + 73 s = 
2/32/42/52/6 + 2/52/62/72/8 + 74S = 
2/22/32/42/52/6 + 2/32/42/52/62/7 + 75 s = 
2/22/32/42/52/62/7 + 2/32/42/52/62/72/8 + 2/62/72/82/22/3 + 76 s = 
2/22/32/42/52/62/72/8 + 2/52/62/72/82/22/3 + 77s = 
2/22/32/42/52/62/72/8 - 1 + 78 s = 

The mixed volume for this system equals 25 and is exact. Of the 25 solutions, eight solutions 
have s = and are thus proper solutions. Two of the eight solutions are real. Denoting I = \J— 1, 
we select the solution 



in v f(y,s) 



(28) 



2/2 = + ^,2/3 = --f,2/4 = -1,2/5 = -1 +-f,2/6 = I,y7 = ^ - ^,Vs = 1 - -f- 
Then we look at the first term of the series 



(29) 
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(30) 



To decide whether the solution is isolated or not, we need to find values for the coefficient of the 
second term in the expansion. Substituting the series in the system in f (y) = and selecting the 
lowest order terms in t leads to an overdetermined linear system in the variables. Solving with 
Maple yields 

1 .1 



Z2 



2' 



Z3 



-1 + J, Z4 = 0, Z 5 = -1, Z 6 = 1 - I, Z 7 



Z8 



1. 



(31) 



Substituting the series in f(y) and we see that the result is 0{t 2 
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7.2 cyclic 12-roots 



Extrapolating on the tropism for cyclic 4-roots, we considered v = (+1,-1, +1,-1, +1,-1, +1,-1, 
+1,-1, +1,-1). For this tropism the first term of the Puiseux series expansion 

X\ = t X2 = t~ l (| — |i\/3) 

23 = —t X4 = t" 1 (— | — |i\/3) 

x 5 =t{-\ + ±i\/3) x 6 = t- 1 (i + Ii^3) 

X7 = — t xg = i _1 (— I + |i\/3) 

X9 = t xio = i -1 (| + gi\/3) 

xn = t (5 - |i\/3) £12 = £ _1 (—3 - 2^^) 

makes the system entirely and exactly equal to zero. Because of the symmetry, we have five other 
solution curves of this type. This is an exact certificate that shows cyclic 12-roots has a curve 
of degree two. Note that 0.866025403784439 is close enoug h to y/3/2 for us to recognize. The 
numerical determination of algebraic numbers in general is done via the integer relation detection 
algorithm of [5], see also [3]. 

The initial root that led to this certificate was one of roots of an initial form system with 
mixed volume — after adding one slack variable to make the system square — equal to 49,816. 
Note that this number is much less than the mixed volume of the original system: 500,352. 



8 Conclusions and Future Directions 

In this paper concepts of numerical algebraic geometry were applied — witness sets and endgames 
- to determine the orientation of tropisms which could lead to certificates for proper algebraic 
curves. 

We sketched how a polyhedral method could pick up all proper algebraic curves defined by 
a polynomial system. To extend this to solution sets of any dimension, say k, one would need 
to consider Puiseux series with k free parameters and look for k dimensional cones of tropisms. 
This generalization leads the development of multiparameter polyhedral homotopies which may 
be of independent interest to numerical analysis. 
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